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Abstract 

Electrons in a multielectron bubble in helium form a spherical, two-dimensional system coupled to 
the ripplons at the bubble surface. The electron-ripplon coupling, known to lead to polaronic effects, 
is shown to give rise also to Cooper pairing. A Bardeen-Cooper-Schrieffer (BCS) Hamiltonian arises 
from the analysis of the electron-ripplon interaction in the bubble, and values of the coupling 
strength are obtained for different bubble configurations. The BCS Hamiltonian on the sphere 
is analysed using the Richardson method. We find that although the typical ripplon energies are 
smaller than the splitting between electronic levels, a redistribution of the electron density over the 
electronic levels is energetically favourable as pairing correlations can be enhanced. The density 
of states of the system with pairing correlations is derived. No gap is present, but the density of 
states reveals a strong step-like increase at the pair-breaking energy. This feature of the density of 
states should enable the unambiguous detection of the proposed state with pairing correlations in 
the bubble, through either capacitance spectroscopy or tunneling experiments, and allow to map 
out the phase diagram of the electronic system in the bubble. 
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I. INTRODUCTION 



Spherical shells of charge carriers appear in a multitude of systems, such as multielectron 
bubbles in liquid helium [j], metal nanoshells coating a non-conducting nanograin 2], carbon 
cages and fullerenes. Although the properties of flat two-dimensional systems have been 
widely studied, revealing new physics, the properties of spherical two-dimensional systems 
are much less well-studied. In this paper, we investigate the possibility and the properties 
of Cooper pairing in the spherical geometry. 

The particular spherical two-dimensional system that we focus on is the multielectron 
bubble in liquid helium. When a flat surface of helium is charged with electrons above a 
critical charge density, an instability occurs with the surface opening to subsume a large 
number of electrons forming a bubble. These multielectron bubbles (MEBs) are typically 
micron-sized cavities inside liquid helium, containing a nanometer thin film of electrons on 
the inner surface of the bubble. The cavity is forced open by the Coulomb repulsion of the 
electrons which is balanced by the surface tension of the helium. The equilibrium shape 
of the bubble is spherical, with a radius R determined by the number of electrons and the 
pressure on the helium. 

The bare single electron states on the surface of the spherical bubble are angular mo- 
mentum eigenstates and have discrete energies, characterized by the angular momentum L 
and with degeneracy 2L + 1. At low temperature there will be a well defined Fermi surface 
located at the highest occupied state. Small-amplitude shape oscillations, including surface 
waves, can be quantized as spherical ripplons. The electrons can interact with these ripplons, 
and we will show that this leads to an attractive effective interaction between electrons. 

This paper has two distinct but interwoven parts. In the first part (section II) we discuss 
in detail how the interactions between electrons and ripplons in the bubble can lead to a 
Cooper pairing scenario. The goal of the first part is to show that a BCS-type Hamiltonian 
provides a plausible description of the electronic system in multielectron bubbles at low tem- 
peratures, and to illustrate the relevant values of the parameters of the model Hamiltonian. 
In the second part (sections III-IV) we investigate the general properties of this Hamilto- 
nian using the Richardson solution for the reduced BCS Hamiltonian describing pairing e.g. 
in nanograins. Both the ground state properties and the density of states of a spherical 
two-dimensional BCS system, are derived and discussed. 
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II. COOPER PAIRS ON A SPHERICAL SURFACE 



First, we investigate the three ingredients of the full Hamiltonian: the electronic part, the 
ripplonic part, and the electron-ripplon coupling. Then, we analyze the effective interaction 
between electrons, resulting from both the Coulomb interactions and the ripplon-mediated 
electron-electron interaction. Finally, we arrive at a BCS-type Hamiltonian for the multi- 
electron bubble. 



A. Electrons and ripplons in the bubble 

The spherical 2D electron system — The well-known Hamiltonian of interacting 
electrons in a flat 2D electron gas (2DEG) in the jellium model can be written in second 
quantization as 

k,cr 

+ J2 Yl Yl Vk+ q , CT 4- q , CT 'Ck>'c kiCT . (i) 

q>0 k,<r k',CT' 

where a c^ a create and destroy an electron with wave number k and spin a, and 

& ,2 e 2 1 

tk = ^ k ' Vk = ^Ak (2) 



where m e is the electron mass, A is the surface of the 2D system, e is the permittivity of 
the medium and e is the electron charge. 

The Hamiltonian of the interacting spherical electronic system has a very similar 
form, provided that one uses spherical harmonics Y Lm (9, 0) instead of plane waves as 
the single-particle basis functions. For this purpose we use the operators t\ m and c^ jm 
that create, resp., annihilate an electron in the angular momentum eigenstate (L,m), i.e. 
ipL,m,{G,(fy = 5 / L,m(#,0)- The Hamiltonian of the interacting spherical two-dimensional elec- 
tron gas (S2DEG) becomes 

L,m,cr 

+ VL ^\L,m)(S)(J,n),a^lL',m')^{J-n),a'^ L ', m 'y^ L ^ m ^ (^) 

J>0,n L,m,a L',m' ,a 



3 



where we use the notation J2l m = SlLo S m =-£ anc ^ 

h 2 L(L + 1) e 2 {-1) L 

eL 2m e R 2 ,VL 2eR2L + V [) 

where R is the radius of the sphere. The typical scale of the kinetic energy e\ = h 2 / (m e R 2 ) 
is listed for several bubble sizes and pressures in table I. The operator c^ L m wj n ) a creates a 
spin a electron in a single particle state resulting from adding the angular momenta (L, m) 
and (J,n). Formally, we have 



W),- 2^ JL, \\ 4 7r(2L' + l) W' U ' L ' U I L ' U ) f5) 
x (J,n;L,m\L' } m'} c£, mV , 

where {J,n; L,m\L' ,m') is the Clebsch-Gordan coefficient for combining the angular mo- 
menta (L, m) and (J, n) into a state of angular momentum (L', m'). The role of the momen- 
tum is now taken by the angular momentum. Indeed, taking L/R — > k for large L links the 
result for the spherical case to that for the flat case. This link was noted previously for the 
structure factor of the spherical 2D electron gas p, and its response to a weak magnetic 
field. 

When the Coulomb energy is small compared to the kinetic energy, the electrons fill up 
a Fermi sea of angular momentum states L (with degeneracy 2L + 1) up to a Fermi level 
L = Lp- The level splitting at the Fermi level, Ae = e^ F — e^ F _i, is given in Table I for 
some typical bubbles. The typical value of the distance between electronic levels is THz (or 
mK). 

Note that the J = term is absent in the Coulomb part of the Hamiltonian (jSJ): this 
term is exactly cancelled by the surface tension energy of the helium as shown in Ref . [5| . In 
Hamiltonian for the flat 2D electron gas in jellium, the q = term is absent because it is 
canceled by a homogeneous positive background introduced in the jellium model. Thus the 
surface tension energy in MEBs takes a role similar to the homogenous positive background 
in jellium. 

Ripplons on the bubble surface - Small-amplitude oscillations of the bubble surface 
can be quantized, leading to the concept of spherical ripplons. The ripplon gas (excluding 
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TABLE I: Several typical length and energy scales for the electron-ripplon system in the bubble 
are given in this table, for bubbles with different numbers of electrons and subjected to different 
pressures. The first row lists the bubble radius in microns. The second row lists the kinetic energy 
scale €\ = ti 2 1 \m c R 2 ) of electrons in the bubble, the energy of angular momentum level L being 
€l = eiL(L + l)/2. The third row gives the electronic level splitting Ae = 6i(Lf + 1), at the 
Fermi level, between the subsequent angular momentum levels of the spherical 2D electron gas. 
The fourth row lists the energy scale for ripplons, hw r = h[a/(pR 3 )] 1 ^ 2 in fiK (1 pK corresponds to 
0.1309 MHz). The fifth row gives the strength of the electric field at the bubble surface, pressing the 
electrons to the helium surface and resulting in electron-ripplon coupling. The sixth row provides 
values for the electron-ripplon coupling constant g in expression @. 







N = 10 4 






N = 10 5 






N = 10 6 




p (Pa) 





10 2 


10 4 





10 2 


10 4 





10 2 


10 4 


R (pm) 


1.062 


.5236 


.1709 


4.937 


1.6930 


.5409 


22.93 


5.393 


1.711 


ei (mK) 


.7836 


3.225 


30.27 


.0363 


.3084 


3.022 


1.68xl0~ 3 


.0304 


.3021 


Ae (mK) 


55.64 


229.0 


2149. 


8.127 


69.10 


673.0 


1.191 


21.53 


213.9 


hjj r (pK) 


10.99 


31.76 


170.3 


1.097 


5.463 


30.25 


.1096 


.9611 


5.378 


E (kV/cm) 


63.80 


262.6 


2465. 


29.54 


251.2 


2461. 


13.70 


247.6 


2459. 


g (mK) 


25.87 


124.1 


1909. 


4.933 


49.70 


803.3 


.9813 


20.61 


338.5 



the breathing mode) is described by the Hamiltonian 

#ripl = ^ ^Lat,nflL,rn- (6) 
L>0,m 

The bare ripplonic frequencies for this system, at a pressure p, are (f| 

where a = 3.6 x 10~ 4 J/m 2 is the surface tension of helium, and p = 145 kg/m 3 is its density. 
In the surface tension dominated regime (pR/cr < 1), ujl = uj r L 3 ^ 2 with uj r = [a / (pR 3 )] 1 ^ 2 . 
Typical ripplon frequencies lie in the MHz-GHz (or pK) range. The ripplon Green's function 
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is defined by 

D(L,m;t) = -i(T[A L>m (t)A+ m (0)) 
with A L>m = a LiTn + al_ m , 

where T is the time ordering operator and A^^ m is a sum of ripplon creation and annihilation 
operators. The unperturbed ripplon propagator (corresponding to a system described by 
.f/ripi above) in the frequency domain is 

2ftjjJL 
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{hjj) 2 — (hull) 2 + if] 
where rj is a positive infinitesimal. 

Electron-ripplon interaction - The electron-ripplon interaction can be written as 

Him = 22MjA J>n ^2 ctL,m)®(J,n),*CL,m,cr, (8) 
J,n L,m,a 

where Mj is the electron-ripplon interaction amplitude and ma c^ L m wj n ) cr^L,m,<7 is the 
(J, n) spherical component of the electron density. 

The interaction between the electrons and the ripplons comes about due to the presence 
of an electric field, generated by the electrons themselves and pressing the electrons against 
the helium surface. This is the electric pressing field E = eN/(2R 2 ), directed radially. 
When a ripplon is present, it moves the electrons in the electric field generated by all other 
electrons and this results in an interaction energy. The interaction energy is the product 
of the displacement caused by the ripplon and the electric field, summed for all electrons, 
similarly as in [jj. Rewriting this interaction energy in second quantization operators, we 
find the interaction Hamiltonian (jHJ with the interaction amplitude 

(J + 1/2) 1/2 

Mj = g [(J+l)(J 2 + J + 2 + 2pR/a)} 1 /^ (9) 
where the coupling constant due to the pressing electric field E (see Table I) is given by 



gV = A=—\—- ( 10 ) 

An additional contribution to the interaction energy between electrons and ripplons can be 
derived as the change in polarization energy of the electron-helium system when the helium 



surface is deformed and the electrons are at rest. This mechanism for coupling was first 
derived by M. Cole for electrons on a flat surface. Following the arguments of Cole, we 
obtain a similar expression for the interaction amplitude but with a different coupling 
constant 

„<p) = h^ £ ~ l e " (11) 
y 8 e+ l4nd 2 R]j 2pRu r ) 1 ; 

where d is the expectation value for the distance between the electron and the helium surface. 

The total electron- ripplon coupling constant is then g = + g( e \ 

A single electron coupled to a bath of ripplons forms a ripplonic polaron j^. In a mul- 

tielectron bubble, the electric field pressing the electrons against the helium surface can be 

much larger than the field achievable on a flat helium surface, so that the ripplonic po- 

larons will be in the strong coupling regime, and can even form a Wigner lattice of ripplonic 

polarons jiol ]. 



B. Effective electron-electron interaction 

Cooper's argument - In this subsection, we follow Cooper's argument for pairing ^l] 
and apply this to the present case of electrons and ripplons in the bubble. The effective 
electron-electron interaction is the sum of the Coulomb interaction between the electrons and 
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we drew a Feynman 



a ripplon-mediated attractive interaction between the electrons. In 
diagram to represent the Coulomb interaction between the electrons. Now, we can add 
another diagram with the same electron propagator lines, but instead of exchanging a virtual 
photon, exchanging a virtual ripplon. This is shown in Fig. ^ 

The two diagrams in Fig. ^ can be represented by a single diagram using an effective 
interaction 

^ m; "> = 2^2lVl + A V)*-(tV + ^ <12) 
Both the Coulomb and the ripplon-exchange interaction are given by a product of two 

vertex factors and one virtual particle propagator. The total electron-electron interaction 
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Ripplon 




FIG. 1: Electrons interact through the exchange of virtual photons or ripplons. The vertex con- 
tributions and propagator for both interactions are shown in this figure. 



Hamiltonian can be written as 

= E E E 

L,m,a L' ,m' ,a' J,n 

Let's study for which regimes V c s(J, n;cu) is attractive, i.e. for which values of (J,n;u) the 
ripplonic part dominates and is attractive. It is clear that small energy transfers make the 
ripplonic part attractive because D^(L,u — > 0) = —2/(ftiu L ). Moreover, in general ripplon 
exchange will indeed occur with u = 0. The reason for this is that the ripplonic energies 
are much smaller than the electronic level spacing, as can be seen from comparing rows 3 
and 4 of table I. For ripplons with L smaller than ~ 10 3 , the absorption or emission of a 
ripplon with angular momentum L cannot change the angular momentum of the electron 
due to energy conservation requirements. The effective interaction at uo = is 

p 2 1 9 M 2 

w^mwr^ (13) 

The attractive interaction dominates strongly at small L, since 2g 2 /(hu r ) 3> e 2 /2eR as can 

be checked for typical bubbles by substituting the values from Table I. It is strongest for 

small L and decreases roughly as L~ 2 : 

2M 2 2g 2 (£ + 1/2) 

V M{ L, m; 0) „ = -J- {L+m l + L H + 2pR/a) ■ (14) 
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dw 
2^ 



V eS (J,n; u)e 



iuit 



The strength of the effective interaction is 2g 2 / (hu a ). Since the g is of the order of mK and 
hjj r of the order of /iK (see Table I), the effective attractive interaction is large as compared 
the relevant energy scales g, hui r , and e\. We are clearly in a strong coupling regime, in 
agreement with the results from [lfj. 

Intralevel pairing - Since the ripplon energies are much smaller than the electronic 
level spacing at the Fermi level, the attractive interaction only takes place between two 
electrons on the same angular momentum level, and these electrons will be scattered into 
final states also on that angular momentum level. An electron in angular momentum state 
|initial) = \L,m), which emits a spherical ripplon in angular momentum state \J, n), finds 
itself in the following superposition of angular momentum states 

L 1 — \L — J\ V 

x (L, m; J, —n\L', m — n) \L', m — n). 

The projection of this final state | final) on the angular momentum level L of the initial state 
is 



f OG [{L, m), (J, -n)] = y (L, 0; J, 0\L, 0) 

x (L, m; J, —n\L, m — n) . 

Thus, the scattering process between two electrons with spin a and o' on the angular mo- 
mentum level L can be described in second quantization as 

L L 

iW= E E Y,^G[(L,m),(J,-n)]f CG [(L,m'),(J,n)} 

m=—Lm'=—L J,n 
X V e g(J,Ti; 0)c^ m _ na C^ tm , +na ,CL,m',a'CL,m,cr 

This interaction Hamiltonian derived for the multielectron bubble is already close to a BCS- 
like interaction Hamiltonian. It involves only electrons on the same angular momentum level 
and couples them with an attractive potential. 

An initial state with a pair characterized by {m,m'} can be scattered into a pair with 
{m — n,m' + n} in various ways: namely by the scattering of a virtual ripplon with J = 
n, n + 1, n + 2, .... provided that J is an even mode to obey the triangle rule of addition 



of angular momenta. All these processes are indistinguishable (initial and final states are 
exactly the same) and their diagrams should be added to get the overall amplitude. Different 
combinations of Clebsch-Gordan coefficients will occur, which can take either a positive or a 
negative sign: the different diagrams can interfere constructively but also destructively. The 
only case where we are sure that the different contributions will interfere constructively, is 
for pairs of electrons with opposite angular momentum m' = —m. The reason for this is that 
(L, m; J, n\L, m + n) has the same sign as (L, — m; J, — n\L, —m — n). We investigated this 
point numerically, and found that the total effective interaction is indeed strongly reduced 
for pairs that do not have opposite angular momentum (m! ^ —m). However, for m! = m 
we find that the effective interaction potential is not reduced and is only weakly dependent 
on m. 

Effective Hamiltonian - From the previous section we know that the interaction be- 
tween the electrons can be written in the form of a BCS interaction Hamiltonian given 
by 

L L 

Hint = - ^ ^ ^ Kn,m',LC^_ m , i c2 m ,. T C L)m;T C L _ m;i , (15) 
L m=—Lm'=—L 

where a — f, J, denotes spin up and spin down and 

x /c G [(L, m), (J, rri - m)]f CG [(L, -m), (J, m - m')]. (16) 

Due to the coefficients (L, 0; J, 0\L, 0) the summation cannot run further than 2L and only 
even values of J contribute. The summation starts from n — \m — m'\, or if this is less than 
2, it starts at J = 2. The J = deformation is not taken into account (it is the radius 
of the bubble), and the J = 1 deformation is a uniform translation which cannot couple 
to the internal degrees of freedom. To proceed, we will introduce an averaged interaction 
amplitude at the Fermi angular momentum level: 

G= (2L F \iy H V ^,l f . (17) 

m,m' 

The interaction amplitude still depends on the angular momentum level L. But, as we 
shall see in the next section, pair correlations occur only in levels close to L F where the 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
L°9 10 [P (Pa)] 

FIG. 2: BCS interaction strength G (in mK) as a function of pressure p and number of electrons 
N in the bubble. 

L-dependence of the interaction amplitude can be neglected. So we use the value of the 
interaction amplitude at Lp also for levels close to Lp. Values for G for various configurations 
are given in Figure El 

With this, we have established that the properties of the ripplon-mediated electron- 
electron interaction lead to a Cooper-type attractive interaction between the electrons, and 
to a BCS-like Hamiltonian 



H e fi — ^LC^ m . a CL, m ;a 
L,m,a 

L 

~ G J2 Yl d t,-m'l d tm';^L, m -A- d L,-m;l- (18) 
L m,m'=—L 

The peculiarity of this pairing Hamiltonian is that the pairing takes place within discrete 
energy levels. This effective Hamiltonian for electrons pairing due to an attractive interac- 
tion brought about by ripplon exchange can be solved by introducing a variational many 
body wave function as in the BCS treatment. However, we chose to apply the Richard- 
son method Q], initially developed in the context of nuclear physics and recently reintro- 
duced [3] to the condensed matter community to describe superconductivity in nanosize 
metallic grains. This method is particularly suitable for finite systems with a discrete level 
structure such as the multielectron bubble. 
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III. RICHARDSON MODEL FOR PAIRING IN A S2DEG 



Having argued that multielectron bubbles are a suitable candidate to observe pairing of 
electrons in a spherical two-dimensional system, we will apply the Richardson model to the 
effective pairing Hamiltonian (|18|) to gain insight in the properties of the paired phase. 

Note that the analysis here can be applied to the general problem of a S2DEG with 
attractive interactions between electrons on the same angular momentum level. The previous 
section provides possible values for the coupling constant G (see Fig. |2J) and for the relevant 
energy scales (see Table I), and a justification for the applicability of the effective pairing 
Hamiltonian (J 18)) to multielectron bubbles specifically. Nevertheless the results derived in 
this section can be used to investigate other systems such as the superconducting properties 
of thin electronic nanoshells, or can be investigated as an academic question regarding 
spherical electronic systems. 



A. Energy levels of interacting electrons 



The Richardson model provides a method of solution for the so-called reduced BCS 
Hamiltonian 

#bcs = Yl e Kko H'^'Akv (19) 

i,cr 

Now consider the Hamiltonian (|18|) that we derived in the previous section, and collect all 
terms that contain operators working on the angular momentum level L: 

Hl = Y e LC[ ma c L:m - a 

L 

~ G 6 t-m'| 6 im';t a WCL,-m;|- (20) 

m,m'=— L 

Note that only electrons within the same angular momentum energy level L interact: only 
mtralevel interactions take place, but no mterlevel interactions. This is due to the fact 
that the relevant ripplon energies are much smaller than the interlevel energy splitting (see 
Table I for typical values in multielectron bubbles). So, the set of electrons with a given 
angular momentum L can be considered as an independent subsystem, described by the 
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Hamiltonian (|20|). The full system is the collection of independent subsystems characterized 
by different L. The full Hamiltonian (fTSj) is just the sum of the Hamiltonians (j2*Uj) describing 
independent subsystems with different L : 

H cS = J2Hl- (21) 

L 

Moreover, each of the Hamiltonians (j2*U)l corresponds to a reduced BCS Hamiltonian ()19|) 
that can be solved with the Richardson method. Indeed, setting Wi : £j = and 

&m;t — CL,m;h b m ;l = CL-m;li (22) 

brings ()20j) into the same form as ()19|) . So, we have a collection of independent systems 
(each characterized by a particular value of L) that can each be solved by the Richardson 
method. As shown by Richardson 12], the exact solution of the reduced BCS Hamiltonian 
for n electron pairs amounts to solving a set of n nonlinear coupled equations. In general, 
the aforementioned set of equations can be solved only by numerical computation. However, 
in the particular case when all the involved single-particle states belong to one an the same 



energy level - as is the case for a spherical multie 
can be easily found analytically (see, e.g., Ref. 



ectron bubble - the energy of electron pairs 



171 ). The result for the energy of electrons 



in the subsystem characterized by angular momentum L can be written down as 

E L,n L , gL ,b L = (2n L + b L )e L - G(n L - g L ) (2L - b L + 2 - n L - g L ) (23) 

The energy levels of the subsystem with angular momentum L is characterized by three 
quantum numbers, n^, and gi- Here is the number of electrons pairs, bj, is the 
number of unpaired electrons, and gi of elementary bosonic pair- hole excitations |l6( in the 
system of pairs. 

To better understand these quantum numbers, consider two bare single-electron states 
\L, m; |) and \L, —m; j). If both states are occupied, this represents an electron pair. The 
electron pair can scatter into another pair of states \L, m'; |) , \L, —m'\ [) under the influence 
of the interaction term in the Hamiltonian (|2(Jj). The number of such pairs with given L is n^, 
and it has to be less than or equal to 2L + 1. Now consider the case where only one state of 
the pair \L, m; |) and \L, —m; j) is occupied. Then we have an unpaired electron that cannot 
participate in the scattering described by the interaction term in (J2*U|) . Moreover, electron 
pairs cannot scatter into the pair of states \L,m; 1),\L, — m\ j) because one of these states 
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is already occupied. The states \L, m; ]),\L, —m; [) are then blocked for scattering of pairs. 
The number of these blocked spin-degenerate bare states equals the number of unpaired 
electrons, and is denoted by b^. The total number of electrons in angular momentum level 
L is then 2np + bp and this has to be less than or equal to 2(2L + 1). The last quantum 
number gp indexes the so-called 'pair-hole excitations' Q|. These are bosonic excitations 
that involve a redistribution of the amplitude for correlated pairs over the 2L + 1 states with 
angular momentum L. Note that gp = corresponds to the ground state for the correlated 
pairs. 

The total energy can be expressed as a sum of the energies for each independent subsys- 
tem: 

oo 

E{n L ,b L ,g L } L=lt2 ,3,... =J ^^^L,n L ,g L ,b L - (24) 
L=0 

The state of the entire system is characterized by a large set of quantum numbers, three 
(riL,bp,gp) for each angular momentum subsystem. 

B. Ground state properties 

How do we characterize the ground state of electrons in a spherical bubble with pairing 
interactions ? In bulk BCS superconductors only electrons in an energy band of the Debye 
energy hiop, around the Fermi level participate in the pairing. However, in the present case, 
the relevant ripplon frequencies are much smaller than the splitting between consecutive L 
levels near the Fermi energy. Therefore one might argue that pairing correlations only take 
place in the one subsystem with L = Lp where Lp is the angular momentum at the Fermi 
level. Indeed, all levels above the Fermi level {L > Lp) are empty, and all levels below the 
Fermi level (L < Lp) are completely filled so that no pairing correlations can be achieved 
by ripplon-mediated scattering. 

Yet this turns out to be wrong. The main difference between pairing in the present case 
and pairing in conventional superconductors is that in conventional BCS superconductors the 
gap A is much smaller than the relevant phonon energy (A -C Fu^d), whereas in the present 
case the interaction energy per electron can be larger than the relevant ripplon energies 
(A twj r ) and even larger than the level splitting (A > tp F — tp F -\). This can be inferred 
by comparing the interaction strengths listed in Table II to the ripplon energies and level 
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splittings reported in Table I. It then becomes energetically advantageous to redistribute the 
electrons amongst the energy levels L. Indeed, by promoting a pair of electrons from a level 
Li below the Fermi level {L\ < Lp) to a level L 2 above the Fermi level (L 2 > Lp), both 
these levels L 1; L 2 can also form pairing correlations. The energy gain by forming pairing 
correlations is larger than the energy needed to promote the electrons form level Li to level 
L 2 . So we obtain the remarkable result that for sufficiently large interaction strength G (as 
we think is the case for MEBs), although only intralevel scattering can take place, still levels 
well below and above the Fermi level can be affected. 

To see this in more detail, let us consider for example a bubble with an even number of 
electrons. The ground state in an even bubble is achieved by gi = (no pair- hole excitations) 
and bi = (no unpaired electrons). In order to find the ground-state configuration rii, 
consider the change in the total energy due to a transfer of an electron pair from the L th 
level to the next higher level. Using Eq. with gt = and hi = 0, we find 



where and n^+i are the number of electron pairs on the levels L and L + l, respectively, 
before the transfer. The first term in the right-hand side of (f23|) represents the energy cost 
in promoting a pair of electrons from level L to level L + l, and the second term corresponds 
to the gain in energy due to pairing correlations. As seen from ()25|). such a transfer reduces 
the total energy if the inequality 



is satisfied. This inequality is not satisfied if til < ^l+i, so that an inverted population of 
the levels obviously never appears in the ground state. Condition (ffilj) is also not fulfilled 
at weak interaction, i.e. for G/t\ < (L + l)/(riL — n^ +1 ). In the case of L 3> 1, the above 
definition of weak interaction can be simplified to G < h 2 /(2m e R 2 ). From Table I and 
Table II, we infer that for MEBs it is likely that G > h 2 / (2m e R 2 ) so that a redistribution 
of electron pairs as described in the previous paragraph is indeed energetically favorable. 

Fig. 01 illustrates the ground-state configuration obtained by minimizing the total energy 
(j2U) with respect to the n^'s (and setting bi = gi = 0). The results are shown as a function 
of G/e± for fixed Lp = 26. The left panel shows the result for Np = 52 electrons on the Fermi 



AE LtL+1 = 2ei(L + 1) - 2G(n L - n L+1 ) 



(25) 



G{n L -n L+1 ) > ei(L + 1) 



(26) 
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FIG. 3: Filled and empty states are indicated with blue resp. green and shown as a function of the 
dimensionless interaction strength G/e\. The typical ripplon energies are smaller than the splitting 
between successive L levels (with degeneracy 2L + 1), so that only intralevel scattering takes splace. 
Nevertheless the electrons are redistributed over different L levels around the Fermi level Lp = 26, 
because the pairing energy is comparable to the level splitting. The left panel shows the result for 
Np = 52 and the right panel for Np = 106 electrons on the Fermi level at G = 0. 

level in the ground state of an MEB with G = 0. This corresponds to roughly half-filling of 
the Fermi level. The right panel shows the result for Np = 106, a closed-shell configuration 
in the G = ground state. Blue (green) color corresponds to filled (empty) states on the 
energy levels. Switching off the interactions, G = 0, we see that all levels below Lp = 26 are 
completely filled (blue color), and all levels above Lp are completely empty (green color). 
For Gjt\ = 2.5, some empty states (green) appear in levels L = 24,25 that are below Lp, 
and some electron pairs appear in levels L = 27,28 that are above Lp. Note that G needs 
to exceed a critical value (of the order of ei/2) for the redistribution of electron pairs to take 
place. Also note that approximately 2G/e\ levels around the Fermi level are affected by the 
redistribution of electron pairs. 

Having obtained, for the ground state, the total energy E S _ B _, we can derive the conden- 
sation energy 

E c = E g . s .(0) - E g . s .(G) - G \N/2] , (27) 
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G/e, 

FIG. 4: The condensation energy E c is shown as a function of the interaction strength G, for 
different numbers of electrons in the Fermi level Lp = 26 at G = 0. As the coupling increases, the 
condensation energy becomes less sensitive to the precise filling of the Fermi level. 

where E gs (G) [E gs (0)} is the ground-state energy in the presence [in the absence] of the 
pairing interaction. The last term in the rhs describes the interaction energy for the un- 
corrected Fermi ground state and \x~\ means the integer part of x. In Fig. |3]we plot the 
condensation energy E as a function of G/e\. First consider the regime of strong inter- 
actions, G ^> ex/2. For this regime, we find that the Np— dependence of the condensation 
energy becomes negligible. The condensation energy rapidly rises with G/e\, approximately 
as E c pa (2LpG) 2 /(3ei). Next, consider the regime of weak interactions (G < ei/2) where, 
as we have seen, the electron redistribution between the energy levels is absent. In this case 
the condensation energy is a linear function of G, and such a linear behavior can indeed be 
seen in Fig. Hjat G < e\/2. In the regime of weak interaction, the condensation energy is 
strongly influenced by Np: it is zero for closed shell configurations and reaches a maximum 
at half filling. For typical MEBs, we find that E /ks is of the order of several Kelvin. 
Decreasing the number of electrons or pressurizing the bubble increases the condensation 
energy. Bigger bubbles have a smaller condensation energy. 
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C. Experimental signatures 



How can the pairing correlations be probed experimentally ? What properties will dis- 
tinguish the state with pair correlations from the normal Fermi sea ? In the preceding 
subsection, we have found the ground state for the system described by the Hamiltonian 
H e fi, expression (jTSjl . and in the preceding section we argued, based on Cooper's argument, 
why this Hamiltonian is suitable to describe MEBs. Nevertheless, more exotic correlated 
or magnetic states not described by the Hamiltonian ([18)1 cannot be completely excluded 
as candidates for the true ground state: only experiment will give a decisive answer as to 
what the realized state in the MEB will be. It is therefore very relevant to discuss accessible 
experimental signatures of the correlated many-body state that we have described. 

In metallic nanograms and nanowires, pairing correlations (leading to superconductivity) 

I I 

were revealed through a measurement of the density of states |l8| . Also here, we propose to 
reveal pairing correlations by probing the density of states. In the case of stabilized MEBs 
[19I ] , the density of states can be measured by spectroscopy. Also tunneling experiments are 
possible by placing an electrode close to the bubble, and reducing the thickness of the layer 
of helium between the bubble and the helium. The density of states that would be revealed 
by these experiments is the subject of the next subsection. 

D. Density of states 

In the previous subsection, we found that the ground state is characterized by a redistri- 
bution of electrons between degenerate bare energy levels, as compared to the ground-state 
electron distribution in the absence of interaction. Here we analyze the effect of the pairing 
interaction on the density of (many-electron) states in MEBs. This density of states can be 
written down in general as 



From expression (J2l|) it is clear that the many-electron energy levels Ei are characterized by 
a set of quantum numbers "i" := {til, bi, 5 i l}l=i,2,3,...- The summation runs over all possible 
sets i of quantum numbers. The degeneracy of the many-electron energy level Ei is denoted 




(28) 



18 



by jw. Consider one particular angular momentum state L with a given {rip, bp, g^}. The 
number of ways to place the unpaired electrons, multiplied by the number of pair-hole 
excitations with given g L , is J nL , gL ,b L (see Refs. jig, |l7[) : 

J »**=*^ X {(^J_C£.-*), ^ (29) 

with Cj, the binomial coefficients. The total degeneracy is the product of the degeneracies 
of the independent systems, 

J® = \[jn L , 9L ,H- (30) 

L 

For graphical representation, it is more convenient to consider instead of T>(E) the quantity 

pE+S/2 

V 5 (E) = / dEV(E), (31) 

JE-5/2 

which gives the number of (many-electron) states in the energy range of width 5 around the 
energy E. 

Fig. El gives an example of the calculated T>$(E). The calculations are performed for 
Np = 54 electrons in the Fermi angular momentum level Lp = 26, so that in the ground 
state at G = the Fermi level is approximately half filled. Together with the whole spectrum 
of excitations (displayed in light green) we also show (in dark blue) the intralevel excitations 
from the ground state. 

At G = 0, as shown in Fig. there is a set of excitations that correspond to different 
interlevel electron transitions between the single-electron bare energy levels with L close to 
Lp (obviously, intralevel transitions have zero energy in this case). At G = 0, all the many- 
electron energy levels in the system under consideration are highly degenerate. This is due 
to the numerous possibilities to distribute electrons over single-electron states belonging to 
a partially filled energy level with a given L. 

At small nonzero G, the many-electron energy levels for each angular momentum level are 
split. Excitations involving pair breaking (bp ^ 0) and pair- hole excitations (gL ^ 0) start 
to affect the density of states. The result, shown in Fig. \^p, is the appearance of energy 
bands in the density of states. The first band corresponds to intralevel excitations from the 
ground state, the other bands involve interlevel excitations from one angular momentum 
state to another. 
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FIG. 5: The calculated density of states for a spherical 2D electron system with pairing, obtained 
within the Richardson model, is shown in these figures as a function of the energy above the ground 
state energy. The different panels correspond to different strengths of the attractive electron- 
electron interaction term in the reduced BCS Hamiltonian. 



At this point it is interesting to investigate the gap in the excitation spectrum for intralevel 
excitations. In superconduc ting nanograins, measuring such a spectroscopic gap signals the 
onset of superconductivity [lg]. First, consider a process where the ground state (with til 
pairs and = unpaired electrons) is transformed into a final state with n' L = — 1 and 
b' L = 2 unpaired electrons. We will refer to such a process as a 'pair-breaking excitation'. 
A pair-breaking 'gap' A^ p_b ^ can then be denned as the energy necessary to create a pair- 
breaking excitation from the ground state: 

A L ^ = #L,»£-1,0,2 - E LtnL}0>0 . (32) 

Using fl23D, 

Af" b) = G(2L + l). (33) 

Similarly, we can calculate the smallest energy needed to create a pair-hole excitation (g^ = 
— > gi = 1) and define a pair-hole excitation gap A^ p_h ^ 

= G(2L + 1). (34) 
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So, in the system under consideration, the smallest energy for a pair-breaking excitation 
equals that for a pair-hole excitation, A^ p_b ^ = A*- P_h \ and we can drop the superscripts 
(p-b) and (p-h). At G ~ ei/2 the intralevel excitation gap A approaches the energy spacing 
ei(L + 1) between single-electron bare energy levels with consecutive values of L. This 
means that, for G ^> ei/2, interlevel excitations exist with energy smaller than the gap A. 
Moreover, as implied by Eq. ()25|) . an increase of G can substantially reduce the energies of 
interlevel transitions of pairs. At G = 0.625ei (see Fig. (St), the lowest interlevel excitations 
(green) already have energies smaller than A (first nonzero blue line). 

At even larger G, pairing correlations appear on many different angular momentum lev- 
els. The pair-breaking/pair- hole excitation gap, A^, depends on the angular momentum. 
Therefore, as seen from Fig. EJi, the peaks of Vs(E) corresponding to intralevel excitations 
from the ground state, are split into peaks corresponding to different angular momentum 
states where pairing can take place. With increasing G, the excitation spectrum tends to 
become quasi-continuous, with jumps of several orders of magnitude in T>$(E) near the en- 
ergies corresponding to pair-breaking or pair- hole excitations. Between these jumps there is 
a more uniform distribution of excitations, corresponding to interlevel transitions of pairs. 

Fig. El provides an "overview" of the behavior of T>s(E) as a function of both E and G for 
fixed Lp = 26 and two different values of Np. One can see an interplay between intralevel 
excitations, whose energies always increase with increasing G, and excitations corresponding 
to interlevel transitions of pairs. For the latter, both an increase and a decrease in energy 
are possible with increasing G. While at G < ei the patterns of T>s(E, G) are very different 
for different Np, at larger G the behavior of T>$(E) in MEBs with a definite parity of the 
number of electrons becomes almost independent of the precise value of this number [cf. 
Figs. Efa) andEfb)]. This "universal behavior" of T>$(E) at large G/ei is further illustrated 
by Fig. El where the calculated dependence of T>$ on (E — E g s )/ ei is shown for different N F 
and Lp. As seen from Fig. El neither moderate changes of Lp nor variations of Np at a fixed 
Lp significantly affect the shape of T>$ versus (E — E g _ B _)/ei for MEBs with a given parity 
of the number of electrons. At the same time, Fig. El demonstrates a pronounced difference 
between the results for even MEBs and those for odd MEBs. 

Summarizing this subsection, we find that in the regime most relevant for multielectron 
bubbles, namely G ^> ei, the pairing correlations reveal themselves in the density of states 
not as a spectroscopic gap, but rather as a significant jump in the density of states at 
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FIG. 6: Density of states of the spherical 2D electron system with pairing interactions, as a function 
of the interaction strength and the energy above the ground state, both expressed in energy units 
€\. The results are shown for two different values of Np, the number of electrons on the Fermi level 
at G = 0: (a) Np = 54 (approximately half filling) and (b) Np = 106 (closed-shell configuration). 

the pair-breaking energy A. Moreover, in this regime {G 3> ei) the density of states is less 
sensitive to even-odd effects. The observation (through tunneling or spectroscopy) of a jump 
in the density of states at A can be used as a way to infer the presence of pairing correlations 
in the MEB. Moreover, the pair-breaking energy can be used to estimate a temperature T c 
above which the pairing correlations will be suppressed: when ksT > A the thermal energy 
is large enough to support an appreciable amount of pair-breaking excitations. For typical 
MEBs, this temperature is of the order of Kelvins. Smaller bubbles or compressed bubbles 
have a larger gap and thus a larger T c . 
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FIG. 7: Density of states as a function of the energy above the ground state energy for MEBs 
with the pairing-interaction strength G = 2e± and different number of electrons. Panel (a) shows 
the results for Lp = 26 and different values of Np , which determines the filling of the Fermi level 
at G = 0. Panel (b) shows the results for two different Lp and different parity of the number of 
electrons in an MEB in the case when the Fermi level is approximately half-filled at G = 0. 

IV. CONCLUSIONS 

In the first part of this paper, we have analyzed how the electron-ripplon interaction on 
a spherical surface may lead to an attractive effective electron-electron interaction and give 
rise to a Cooper pairing scenario. The effective Hamiltonian of the two dimensional spherical 
electron system is mapped on a BCS-type Hamiltonian and typical values of energies, length 
scales and interaction strengths are estimated. 

In the second part we use Richardson's method to investigate pairing properties of a 
two dimensional spherical electron system. We find that when the condensation energy per 
pair is larger than the bubble energy scale h 2 /(m e R 2 ), the ground state of the system ac- 
quires unique properties that set it apart from pairing in conventional superconductors or 
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superconducting nanograms. In particular, we show that although only intralevel interac- 
tions are included (since the relevant ripplon energies are smaller than the level splitting), 
electron pairs nevertheless redistribute themselves among the different levels and pairing 
takes place in an interval of energies around the Fermi energy, much larger than the typ- 
ical ripplon energy. The density of states reveals an intricate interplay between intralevel 
transitions and interlevel excitation of pairs, evolving from a discrete spectrum typical for 
confined systems to a staircase quasi-continuum upon increasing interaction strength. At 
strong coupling the density of states reveal the presence of pairing correlations not through 
a spectroscopic gap as in metallic nanograins 3], but through a jump in the density of 
states at the pair-breaking energy. 

These results show that spherical electron systems reveal particularly interesting pairing 
properties, distinct from their bulk or flat-surface counterparts, combining both topologi- 
cal effects and confinement effects. MEBs are a particularly pure realization of the two- 
dimensional electron system (just as an electron film on helium forms a pure realization of a 
flat 2DEG). Thin nanoshells (monolayer gold coatings of a non-conducting nanograin) may 
be another realization of this system. 
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